Radiative decays of charmonium into light mesons 
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We apply perturbative QCD to the radiative decays of charmonia J/tp and Xcj into light mesons. 
We perform a complete numerical calculation for the quark-gluon loop diagrams involved in these 
processes. The calculated J/tp decay branching ratios into P-wave mesons /2(1270) and /i(1285) 
fit the data well, while that of /o(980) (if treated as an ss meson) is predicted to be 1.6 x 10 -4 , 
which implies that fo (1710) can not be the ss or (mm + dd) /\/2 meson. Decays of P-wave charmonia 
XcJ — * p( w ,0)7 (J=0,l,2) are also studied, and the branching ratio of \a — > P7 is predicted to be 
1.4 x 10~ 5 , which may be tested by CLEO-c and BESIII with future experiments. 

Decays of heavy quarkonium into light hadrons are forbidden processes by the OZI rules. The radiative decays of 
J/tp into light hadrons such as J/tp — > 7/2(1270) are expected to proceed via two virtual gluons that subsequently 
convert to light mesons, with the photon emitted by the charm quarks (see Fig.l for the three Feynman diagrams). 
These processes are interesting since they can provide useful information for understanding the fundamental theory 
of strong interactions, namely, quantum chromodynamics (QCD). 

In this letter, we will report calculations for radiative decays of charmonia J/tp and Xcj into various light mesons 
based on nonrelativistic quantum chromodynamics. Hopefully, the J/tp mass is large enough for using perturbative 
QCD to calculate these processes. 

Differing from previous works (see, e.g.,[l|,[2|), we will perform a complete numerical calculation for the quark-gluon 
qq \ loop diagrams involved in these processes. 

. Here, it is emphasized that in our approach the description of J/tp as a nonrelativistic cc bound state is a good 
approximation; while treating the light mesons such as /2(1270) as nonrelativistic qq bound states (q = u,d,s) with 
constituent quark masses m u — vrid ~ 350 MeV, m s ss 500 MeV is only an approximation for describing the 
nonperturbative dynamics that has been used in the nonrelativistic constituent quark models rather successfully, and 
was also used in previous studies for these radiative decay processes 0, 0]. In fact, if the current quark masses for the 
light quarks, which will vanish in the chiral limit, were used in the quark-gluon loop diagram, we would encounter 
difficulties with soft and collinear divergences. 
£pH, In the following, we will first focus on the radiative decays of J/tp to P-wave light mesons such as /2(1270), /i(1285), 
O ' and /o(980). We adopt the descriptions of nonrelativistic bound states, i.e., both J/tp and light mesons are described 
by the color-singlet non-relativistic wave functions. As widely accepted assignments we assume that ^2(1270) and 
' /i(1285) are mainly composed of (uu + dd)/\f 7 i (neglecting the mixing with ss for simplicity). But for /o(980), there 
, are many possible assignments such as the tetraquark state, the K K molecule, and the P-wave ss dominated state 
J^~j ■ (for related discussions on Jo (980) and other scalar mesons, see, e.g., the topical review-note on scalar mesons in Q 
and [H). Since experimental data show that Df — ► fa^ + has a large branching ratio (BR)[3[, we prefer to assign 
■ Jo (980) as an ss dominated P-wave state as a tentative choice. 

In the calculations, the spin projection operators for 3 Si, 1 S$ and 3 Pj states and their relations can be found in 
[Eil- We use FeynCalc 7] for the tensor reduction and LoopTools[8] for the numerical evaluation of infrared safe 
integrals. We follow the way in Ref. @ to deal with five-point functions and high tensor loop integrals that can not 
be calculated by LoopTools and FeynCalc, such as 

Eaf}pa = I dDk ^[(k+p q y - m 2](fe + 2p q Y[{p c + fc) 2 - m2][(fc + 2p q -p c )2 - ml] ' (1) 

where p c is the momentum of c quark, p q the momentum of light quark, m c the charm quark mass, and m q the 
constituent light quark mass of fj meson. For the constituent quark masses, we choose m u — rrid — 0.35 GcV, 
m s = 0.50 Gev, m c = 1.5 GeV. 

Using the theoretical expression for the widths of /a ^77, 

r (th) 6JV e a 2^ n . 2 l^p(0)l 2 r («o _12N i 2 \U' P (0)\ 2 m 

where N c =3 is the color number, a=j|^, and fitting them with their experimental values 2.6 KeV and 0.081 KeV for 
/ 2 (1270) and /£(1525) respectively 0, we get 

l^/ 2 (i270)(0)| 2 =6.6xl0- 4 GeV 5 , |^, (1525) (0)| 2 = 1.1 x lO" 3 GeV 5 . (3) 
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FIG. 1: Feynman diagrams for charmonium ( 3 <Si) radiative decays into light mesons ( 3 Pj) 
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TABLE I: Numerical results for B{J/ip -> /,/7). 



For charmonium, the value of the wave function at the origin is taken to be [ioj |: 

\K J/4 ,{0)\ 2 = 0.81 GeF 3 , |^ c/ (0)| 2 = 0.075 Ge^ 5 . (4) 

In the calculation, the meson masses are taken to be M = 2m q in the nonrelativistic limit. For the J ftp decay we 
take a s (2m c ) = 0.26. Our results are listed in the Table I. 

Then we will give the contributions of different helicities. As in Ref.[ll|, we choose the moving direction of fj as 
the z-axis, and introduce three polarization vectors: 

= ^=(0,1,1,0), ^(-l)~(Q,l,-i,0), ^(0) = l(|k|,Q,0,fc°), (5) 

thus we can characterize the polarization tensor e tlv (\) of /V With the same parameters, we give the branching ratios 
for different helicity states in Table II. The corresponding values of helicity parameters x and y are x = 0.46, y = 0.23. 
Recently, new experimental data for the contributions of different helicities in process J ftp — > 7/2(1270) have been 
given by the BES Collaboration [12]: x = 0.89 ± 0.02 ± 0.10 and y = 0.46 ± 0.02 ± 0.17 (see also 01). It is about 
2 times larger than our results. But if we use a larger constituent quark mass, e.g. m u — M(1270)/2, we will get 
substantially increased values x = 0.79 and y — 0.58 (also see Ref.[l(). We emphasize that the helicity parameters 
are very sensitive to the light quark masses, and hence very useful in clarifying the decay mechanisms. Note that if 
Tin/wc — > 0, we will have x — > and y — > 0, but this is inconsistent with data. 

From Table I we sec BR of J/t/j — > /o(980)7 (assuming /g(980) to be ss) is predicted to be 1.6 x 10~ 4 . Actually, this 
calculation can also apply for /o(1710) (if treated as ss) with an even smaller BR due to the phase space correction. 
However, experimentally BR of J/ip — > /q(1710)7 is larger than 9 x 10~ 4 [3]. This implies that the /o(1710) can not 
be a pure ss or (uu + dd)/\/2 scalar meson (for the latter the BR is predicted to be about 4 x 10~ 4 ). One possible 
explanation is that the /o(1710) has a substantial gluonic component due to qq-ghieball mixing. 

Next, we discuss J/ip radiative decays to 77 and rf . Assume rj' and r\ can be written as admixtures of flavor-singlet 
ipi and octet ip$ : 

ipi = —=(uu + dd + ss), -08 = —={uu + dd — 2ss), (6) 
V3 V6 

7] = -08 cos 6 — ipi sin 0, 77' = -08 sin + 01 cos 8, (7) 
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TABLE II: Results for J/0 — » /,/7 with different helicity states. 
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TABLE III: Branching ratios for J/ip — > T](r)')~/. 
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TABLE IV: Predicted branching ratios for XcJ — * r yp(&, ( />)- 



where 9 is the mixing angel, taken to be a widely accepted value 6 = — 20°[lJ|. For the 77(77')) the wave function at the 
origin may be determined from the leptonic decay of vector mesons V — ► e + e~ (V = tf>, p), with the assumption that 
in the nonrelativistic quark model the wave functions of the *~ and 1 mesons should be equal in the nonrelativistic 
limit. Using 



T(0(1O2O) -> e+e") = N c Q 2 s a 2 = (1.27 ± 0.04) KeV, (8) 

we can get |7?.(0)| 2 =0.054 GeV 3 for ss. Similarly we can get |7?.(0)| 2 =0.032 GeV 3 for uu and dd from p leptonic decay. 
The decay branching ratios of J/ip — ► 777(77') arc then calculated and listed in Tabl dllll 

We can see that for J/tp — ► 777(77') the simple nonrelativistic quark model calculations fit the data rather well, 
and the values are compatible with other approaches e.g. the 77(77') — r\ c mixing mechanism based on the QCD axial 
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Another interesting processes are decays of P-wave charmonia \cJ ~ * p( tL '?'/ , )7 (J = Q>1)2), where the photon is 
emitted by the light quarks but not the charm quarks, in contrast to the above discussed J/tp radiative decays. The 
calculated branching ratios are shown in Table HVl In particular, the branching ratio of \ci — * PI is predicted to be 
1.4 x 10~ 5 , which may be tested by CLEO-c and BESIII in future experiments. 

In this letter, we apply perturbative QCD to the radiative decays of J/ip and \cj into light mesons by treating 
the light mesons such as /2(1270) as nonrelativistic qq states as in the nonrelativistic constituent quark models. We 
perform a complete numerical calculation for the quark-gluon loop diagrams, where the light quark masses are taken 
to be the constituent quark masses. The calculated J/ip decay branching ratios to the P-wave mesons /2(1270) and 
/i(1285) fit the data well (but the observed helicity parameter x for /2(1270) favors even larger light quark masses). 
Moreover, the results for S-wave mesons 77(77') seem to be also rather good. Our calculation does not support /o(1710) 
being the ss or (uu + dd)/\/2 meson, while the possible ss interpretation for /o(980) can be tested. The predicted 
branching ratios for XcJ p{^,4>)l (J = 0>1,2) are generally small, and may be tested by CLEO-c and BESIII with 
future experiments. 

However, we must bear in mind that for J/ip decays the applicability of perturbative QCD is marginal, and the 
nonperturbative effects should be important. Moreover, the successes of nonrelativistic constituent quark models for 
light hadrons have not been fully understood from QCD. The constituent quark model description for 77 and 77' is 
even more questionable. In this paper, to evaluate the quark-gluon loop diagrams and to avoid the soft and collinear 
divergences caused by massless light quarks, we have tried to use the constituent quark masses as in the nonrelativistic 
quark model for light mesons. Hopefully, the good (if not accidental) fit of the calculation presented here with the 
experimental data can shed light on a deeper understanding of J/ip decays in QCD, which evidently deserve further 
studies. 
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